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Abstract 

We present some new mathematical tools which help to derive information about the quark 
mass matrices directly from experimental data and to elucidate the structure of these mass 
matrices. 
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One hopes that a successful choice of a mass matrix might help to suggest an ap- 
propriate dynamical model explaining the quark mass spectrum. In any case, one has 
to be able to write a mass matrix consistent with the experiment, regardless of what 
particular dynamics are responsible for generating quark masses. Over the years, a 
number of proposals have been advanced for quark mass matrices, such as those in refs. 

||, |3|, [|, [|, [| . In this paper we would like to present some useful mathematical tools 
for the study of quark mass matrices and the connection with quark mixing. The start- 
ing point of our study is an observation which has so far not been exploited: since the 
Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix V is an element of the Lie 
group U(3), it is in one-to-one correspondence with an element H of the Lie algebra 
u(3). H is thus the infinitesimal generator for V . It is of interest to determine the H 
corresponding to the empirically measured V. As we shall show, one can also use H to 
facilitate the determination of families of quark mass matrices which are consistent with 
experiment. (Here we use the fact [7j that there are three light neutrinos, hence three 
usual SU(2) doublets of leptons, and, as required, e.g., by anomaly cancellation, thus 
also three SU(2) doublets of quarks, so that V G U(3).) 

We denote the quark mass matrices for the up and down sectors M u and corre- 
spondingly. By means of a bilinear transformation one can diagonalize M u and Md- 



where U u ,l-, Ud,L, U u> r, Ud,R are unitary matrices. Below we shall suppress the subscript 
L. 

In the standard model with only left-handed charged weak currents, the matrix Ur 
is not a measurable quantity; and therefore one may choose the U U) r and U^r in such 
a way that M u and Md are both hermitian. Indeed, an arbitrary complex 3x3 matrix 
is defined by 18 real parameters. In order to make it hermitian, one has to impose 9 
constraints to eliminate 9 degrees of freedom. To do that one may use 9 degrees of 
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freedom associated with each of the matrices U u% r and U d ,R. Below, we assume both M u 
and M d are hermitian. 

The experimental data |7j provides one with the absolute values of the matrix el- 
ements of the CKM mixing matrix V = U u U d (taken to be 3 x 3). One also obtains 
certain phase information from experiment which, for a given parametrization, can be 
translated into a determination of the phases of the elements of V. Suppose the pair of 
hermitian mass matrices M u and M d , 



UlD u U u 




(2) 
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U\D d U d 



corresponds to some particular matrix V. Then one can see that for every unitary matrix 
Uq there exists another pair of hermitian matrices: 



' M' u = (U u U o y D u (U u Uo) 

(3) 

k M' d = (U d U o y D d (U d U ) 

that corresponds to the same matrix V. One can also see that the family of matrices M' u 
and M' d as a function of an arbitrary unitary matrix Uq exhausts all possible solutions 
for the mass matrices. This proves that for every mixing matrix V G U(3) there exists 
a 9-parameter family of matrices that might be chosen to be mass matrices. 

Empirically, the data constrains \ V\, the matrix of magnitudes of each of the elements 
of V, to be close to the identity matrix. Therefore, it is useful to represent V as 

V = e iaH = U U U\ (4) 

where H is some hermitian matrix and a is a real number. One may choose H to have 
its dominant (largest, in absolute value) eigenvalue to be of the order of 1, yet smaller 
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than 1. Then for a consistent with the data |7j, we find that |a| ~ 0.3. We will use the 
fact that a < 1 below. 

Suppose that V is known. Six parameters, the phases of the quark fields, are ar- 
bitrary. We use three of them to rephase V to make it be close to the 3x3 identity 
matrix. This is possible since \V\ ~ 1. Then one can calculate the matrix aH using the 
Sylvester's theorem ||: 

The latter is valid as long as all the eigenvalues of V are distinct. If the matrix V has 2 
or 3 degenerate eigenvalues, then ln(V) is not uniquely defined, and there is a family of 
matrices H that satisfy (|j). This family of solutions is explicitly shown in ||. However, 
this mathematical subtlety is not of physical importance since one can always avoid the 
degenerate case by some small change in the phases of quark fields that removes any 
initial degeneracy in eigenvalues. The theorem ([5]) is particularly powerful since it gives 
an exact result in a closed form. 

We will demonstrate how one could construct a set of mass matrices starting from 
the mixing matrix V . In particular, we choose to describe those mass matrices that can 
be diagonalized by a unitary transformation U u -d generated by an element in the u(3) 
algebra that commutes with H . 

In terms of the matrix H, 



V = 1 + iaH - -a 2 H 2 + ... + —UaH) n + ... (6) 
2 n\ 

One can also decompose an arbitrary U\ and U2, that commute with H, in a similar 
series whose coefficients are proportional to certain real numbers: {x^} and {x^} 
where n = 1, 2, 3... 
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' U u = l + iax^H - \x ( ^a 2 E 2 + ... + ±xM(iaH) n + ... 

(7) 

^ U d = 1 + iax[ d) H - \xfa 2 H 2 + ... + ±x^(iaH) n + ... 
This series is convergent if 



sup \x n \ < oo 

n=l,2,3... 

The constraints = V, = 1 and U d U d = 1 are satisfied simultaneously up to 

the third order in a if the coefficients are chosen as follows: 



f U u = l + iaHx - ^H 2 x 2 + ... 

(8) 

{ U d = l + iaH{x - 1) - ^H 2 [x - l) 2 + ... 

where x is some real parameter. Then, to the second order in a, one may express M u 
and Md as functions of x: 

' M u = W U D U U U = 
D u + iax[D u , H] - \a 2 x 2 [[D u , H],H} + ... 

(9) 

M d = U\D d U d = 

D d + ia(x - l)[D d ; H] - \a 2 {x - l) 2 [[D d , H], H] + ... 

We see that, in our parametrization, one gets a good approximation of M's as a 
function of one real parameter. Thus, from V and the resultant H, we have obtained a 
family of (experimentally acceptable) mass matrices depending on a parameter x. 

In passing, we note the following technical point: for rapid convergence of the series 
(9), V should be close to 1, the identity. The three phases of the diagonal elements of V 
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are not physical and they may be changed by rephasing the quark fields. Such a change 
in the overall phase of V corresponds to shifting a by a real number. Clearly, |a| is the 
smallest for V ~ 1, which yields rapid convergence of both (6) and (9). 

The family of solutions (9) for the "up" and "down" quark mass matrices (accurate 
to second order in a) is 

M u « D u + tax[D u , H] - \a 2 x 2 [[D u , H], H] 

(10) 

M d « D d + ia{x - l)[D d , H] - \a\x - l) 2 [[D d , H], H] 
To illustrate the method with a simple numerical example let us take 



V 



/ 0.97525 0.22104 0.00506146"* 40 00 ° 

-0.22099e i0 0359 ° 0.974306"^ 000184 ° 0.043619 
V 0.00666576"* 28 40 ° -0.043403^ a9485 ° 0.99904 



cb\ 




(consistent with current data 0). Note that this gives \V u b\ 
parametrization-invariant CP violation quantity jnj J = /^{VnV^V^V^} the value 
J = 3.1 x 10" 5 . 

First, we calculate the eigenvalues of V: v x = 0.97415 + 0.22587i, v 2 = 0.97443 - 
0.22468z, v 3 = 0.999999 + 0.001221 Next, we substitute these into © to calculate H 
and obtain: 



aH 



I -0.223? 

0.223i 
V -0.0032 + 0.001* 0.044i 



-0.0032 - 0.001* 
-0.044i 




(12) 



We now choose the quark masses to be equal to their running masses evaluated at 
Qo = 1 GeV and with A = 200 MeV in the MS renormalization scheme. From the 
analysis of §: m u (Q ) = (5.1 ± 1.5) MeV, rh d (Q ) = (8.9 ± 2.6) MeV, rh s (Q ) = 
(175 ± 55) MeV, m c (Q ) = (1.35 ± 0.05) GeV, m b (Q ) = (5.9 ± 0.1) GeV. 
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We also use the value m t (Qo) = 220 GeV corresponding to a physical mass of mf ys = 
130 GeV, consistent with current constraints. Specifically, we choose the diagonal mass 
matrices to be 



D u = diag(5 MeV, -1.35 GeV, 220 GeV) , . 

D d = diag(8.9 MeV, -175 MeV, 5.9 GeV) ^ ' 

Suppose x = 1/2. In this case, one retains a symmetry between the "up" and "down" 
sectors in (15). For these values of quark masses the first order approximation of the 
"up" and "down" sector mass matrices gives the following: 



M u = 220GeV 



I -5.09 x 10~ 5 (6.77 x 10- 4 )e- J 3 ° (2.37 x lO" 3 )^ 42 ° \ 

(6.77 x 10~ 4 y 3 ° -5.58 x 10~ 3 -0.0221 
V (2.37 x 10- 3 )e" i 42 ° -0.0221 0.99951 



(14) 



M d = 5.9GeV 



I 1.49 x 10~ 3 -3.47 x 10- 3 (1.78 x lO- 3 )^ 66 ° \ 

-3.47 x 10~ 3 -0.0293 0.0226 

V (1.78 x 10- 3 y 66 ° 0.0226 0.99951 J 



(15) 



Even though this is only a second order approximation it still shows what is the 
relative value of different terms in the mass matrices for this particular case of x = (1/2). 
For instance, the matrix V calculated back from M u and M d is different from the original 
V only in the third significant digit. 

It is of interest to express the CP violation quantity J directly in terms of H. Recall 
that 

J = -^ W det(C) (16) 

where 

F u = (m t - m c )(m t - m u )(m c - m u ) . , 

Fd = (rn b - m s )(m b - m d )(m s - m d ) 
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Now 



C = -i[M U) M d ] 



(18) 



f M u = UlD u U u = W d VW u VU c 



d — 



U\(D U + ia[D u , H] - %[[D W H],H] + ...)U d 



(19) 



so that up to the first non- vanishing order in a, 



iC = [M u , M d ] = i(a + 0(a 2 )) U\[[H, D u ], D d ]U d 



(20) 



Using (pOD , we thus find 



J= (a + O^ 2 )) 3 ^^^^!] 



(21) 



This result allows one to calculate the parameter J up to the small corrections of 
the order of a 2 . The higher order terms can be obtained by retaining the corresponding 



yields J ~ 3.2 x 10~ 5 , very close to the exact value J = 3.1 x 10~ 5 . 

These results thus establish a useful connection between the M matrices, the uni- 
tary transformations which diagonalize them, and the actual quark mixing matrix V. 
Furthermore, they introduce an interesting object, the matrix H in Lie algebra u(3) 
corresponding to V in the Lie group U(3). We have shown how with our tools, one can 
obtain a family of experimentally acceptable mass matrices directly, given knowledge of 
V. This technique complements the earlier typical one of hypothesizing some ansatz for 
a mass matrix and then checking to see if it fits experiment. 



higher order terms in (p0|). For exampli 



.e, for the numerical illustration (11)-(15), (pl|) 
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